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Abstract 

Modified Altarelli-Parisi splitting functions were recently proposed to model 
multi-parton radiation in a dense medium and describe jet quenching, one 
of the most striking features of heavy-ion collisions. We implement medium- 
modified splitting functions in the HERWIG parton shower algorithm, 
which satisfies the angular ordering prescription, and present a few parton- 
level results, such as transverse momentum, angle and energy-fraction 
distributions, which exhibit remarkable medium-induced effects. We also 
comment on the comparison with respect to the results yielded by other 
implementations of medium-modified splitting functions in the framework of 
virtuality-ordered parton cascades. 



1 Introduction 



Measurements performed at the Relativistic Heavy Ion Collider (RHIC) have empha- 
sized the jet-quenching phenomenon, namely the suppression of particle production at 
large transverse momentum (pt) with respect to the naive vacuum expectations, as well 
as other related phenomena, such as the disappearance or distortion of the particles 
directed in opposite azimuth to a triggered high-transverse-momentum one [1]. A typ- 
ical explanation of jet quenching consists in assuming a higher radiative energy loss in 
a dense medium, which, with respect to the vacuum case, allows partons potentially 
produced at large pr to further emit, in such a way to decrease the high-p^ multiplicity 
and enhance the low-p^ one. A lot of work has been undertaken in order to theoretically 
describe such effects [2]. Calculations carried out so far are suitable to predict inclusive 
quantities, but they are not adequate to describe exclusive final states, which is com- 
pelling for the experimental studies. Moreover, they deal with multi-parton emissions 
only by means of simple assumptions [3,4]. 

As happens when studying hadron collisions in the vacuum, Monte Carlo event gen- 
erators are the best possible tool for the sake of performing experimental analyses. In 
fact, Monte Carlo codes, such as the multi-purpose HERWIG [5] and PYTHIA [6] gen- 
erators, contain fairly large libraries of hard-scattering processes, as well as initial- and 
final-state parton showers, models for hadronization and underlying event. Moreover, 
these programs provide one with exclusive final states: the user can set acceptance 
cuts on final-state particles and interface the Monte Carlo output with detector simu- 
lation programs. For such reasons, having Monte Carlo programs capable of simulating 
nucleus- nucleus interactions will be of great interest for the heavy-ion community work- 
ing at RHIC [1] and ultimately at the LHC (see Refs. [7-9] from ALICE, CMS and 
ATLAS experiments, respectively). 

Although an algorithm for multiple radiation in a dense medium should require sev- 
eral a priori assumptions, it was recently proposed [10, 11] that a simple prescription 
to implement medium-induced effects in parton shower simulations consists in adding 
to the Altarelli-Parisi splitting function a term depending on the parameters which 
characterize the medium, on the virtuality and on the energy of the branching parton. 
Under reasonable hypotheses on the medium properties, it was shown in [10] that, even 
allowing only one medium-modified splitting, it was possible to obtain azimuthal dis- 
tributions qualitatively similar to the RHIC observations. Modified splitting functions 
have been recently implemented in the framework of the PYTHIA event generator [6]. 
This implementation is known under the name of Q-PYTHIA: the inclusion of modi- 
fied splitting functions is detailed in [12], whereas the Q-PYTHIA fortran code can be 
downloaded from [13]. Although a thorough comparison with RHIC data has not yet 
been performed, the results of Ref. [12] look in qualitative agreement with the expected 
features of heavy-ion collisions, namely suppression of particle production at large pt, 
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broader angular distributions and enhancement of intra-jet multiplicities. 

In this paper we extend the method discussed in [10-12] and apply it to the HER- 
WIG generator [5], the other general-purpose Monte Carlo program pretty much used 
by the experimental collaborations. In fact, having at least two codes implementing 
medium-modified splitting functions is mandatory for the sake of comparison and esti- 
mating the Monte Carlo uncertainty on a prediction. Also, it is well known that the 
algorithms implemented in HERWIG and PYTHIA are indeed quite different, in both 
parton showers and hadronization, which makes compelling including medium effects 
even in HERWIG and using it to analyse heavy-ion collisions. 

In fact, HERWIG showers satisfy angular ordering of multiple soft radiation, which 
correctly accounts for colour coherence and allows one to probabilistically implement 
multiple soft emissions [14]. On the contrary, PYTHIA cascades are traditionally or- 
dered according to the virtuality of the splitting parton, with an option to reject non- 
angular-ordered emissions. This was the scenario in which the authors of Ref. [12] 
implemented modified splitting functions. The PYTHIA evolution variable is not, how- 
ever, entirely equivalent to angular ordering in the soft limit: although in several cases 
the actual ordering variable does not make really big changes, when comparing with 
experimental observables sensitive to colour coherence, as done in Ref. [15], HERWIG 
agrees with the data better than PYTHIA. Transverse-momentum ordering, included 
in the latest PYTHIA version [16], should yield a better description of angular order- 
ing. Nevertheless, as discussed in [17], discrepancies with respect to HERWIG are still 
present when considering, e.g., the so-called non-global observables [18], sensitive to ra- 
diation in a limited part of phase space, such as the transverse energy flow in a rapidity 
gap. Even for the purpose of hadronization, the two programs implement very different 
models, namely the cluster model [19] (HERWIG), based on colour preconfinement, and 
the string model [20] (PYTHIA), both depending on a few parameters to be fitted to 
experimental data. It is only after turning hadronization on and fitting such models 
to the same data set (see, e.g. the study [21] for heavy-quark fragmentation) that the 
comparison between HERWIG and PYTHIA can be consistently made. At parton-level, 
due to the above mentioned differences in the treatment of parton showers, discrepan- 
cies between HERWIG and PYTHIA must instead be expected, in both vacuum and 
medium- modified cascades. 

More generally, it is worthwhile pointing out that the choice of the ordering variable 
for multi-parton radiation in a medium is currently an open issue (see the discussion 
in [11]). For single-gluon emission, it is known that finite medium-size effects affect the 
radiation pattern by producing a suppression of small-angle emission, equivalent to the 
Landau- Pomeranchuk-Migdal effect [22]. Similar mechanisms may, for example, modify 
the suppression associated with colour coherence, which leads to angular ordering. Im- 
plementing medium-modified splitting functions in angular-ordered parton showers, as 
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we shall do in this paper, or in virtuality-ordered cascades, as done in [12], will allow one 
to study how the jet structure in a medium depends on the adoption of a given evolution 
variable. Other approaches have also been pursued, following analytic or Monte Carlo 
methods [23] ; the experience gained by applying all these different calculations and algo- 
rithms will be essential for a correct characterization of the underlying dynamics when 
experimental data from heavy-ion collisions will become available. 

Hereafter, we shall employ the latest fortran version of HERWIG, with a few routines 
modified to include medium effects. The documentation of this code is currently in 
progress [24]: in the near future, the heavy-ion community will be able to compare 
the results of Q-PYTHIA with the ones of the medium-modified Q-HERWIG code. 
Furthermore, our work can be straightforwardly implemented even in the object-oriented 
C++ versions of HERWIG [25] and PYTHIA [26], containing the basic physics of the 
fortran ones, plus a number of remarkable improvements. 

The plan of the present paper is the following. In section 2 we shall shortly review the 
basics of the HERWIG showering algorithm in the vacuum and how it differs with respect 
to the PYTHIA one. In section 3 we will discuss the implementation of the modified 
Altarelli-Parisi splitting functions and present results for the HERWIG Sudakov form 
factor. In section 4 we shall present a few parton-level results, showing the role played 
by medium effects. In section 5 we will summarize the main points of our investigation 
and make remarks on possible extensions of the analysis here presented. 



2 The HERWIG parton shower algorithm 

In this section we discuss the HERWIG algorithm in the vacuum, to which we shall apply 
later on the modifications due to a dense medium. Monte Carlo algorithms rely on the 
factorization of the branching probability for soft or coUinear parton radiation. Referring 
to final-state radiation, the one mostly affected by medium effects, the probability of 
emission of a soft/collinear parton reads [14]: 



27r ' ' AsiQ^Q 



In ([T]), P{z) is the Altarelli-Parisi splitting function, z the energy fraction of the radiated 
parton with respect to the emitter and is the ordering variable of the shower, whose 
maximum and minimum values are Q^ax Qoy respectively. Q^ax is set by the 
hard-scattering process, whereas Ql is the user-defined scale at which the perturbative 
cascade ends and cluster hadronization starts [19]. 

In HERWIG [5] the evolution variable is = E'^(, where E is the splitting-parton 
energy and C = (pi ■ P2) I {EiE^) , Pi{Ei) and P2{E2) being the momenta (energies) of 
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the radiated partons. For massless partons, turns out to be an energy-weighted 
angle, i.e. ~ E'^iX — cos6), where 6 is the emission angle. For soft radiation, the 
HERWIG evolution variable is thus equivalent to angular ordering [14], valid in soft 
approximation, for azimuthally-averaged quantities, in the large- A'^c limit. The scale of 
the strong coupling constant 05 is the transverse momentum of the radiated parton with 
respect to the emitter. In this way, one includes in the algorithm a class of subleading 
soft /coUinear logarithms, thus improving the accuracy of the shower even beyond the 
leading- logarithmic approximation [27]. 

In ([T]) Asi^Ql, Ql) is the Sudakov form factor, expressing the probability of evolution 
from Ql to Q2 with no resolvable emission; in diagrammatic terms, the Sudakov form 
factor sums up virtual and unresolved real emissions to all orders. The ratio of form 
factors in Eq. ([1]) represents the probability that the emission at {z, Q"^) is the first, i.e. 
that there is no radiation during the evolution between Q^^x Q^- given by the 
following expression: 

As QLx, = exp - / / dz P{z) . (2) 

In Eq. ([2]), the limits of the 2;-integration are given in HERWIG byS z^in = Q^/Q and 

2^max 1 ^min' 

The showering algorithm ([T]) is frame-dependent, however one can prove [14] that, if 
i and j are colour-connected partons, the upper values of the evolution variable satisfy 
the relation (5i,maxQj,max = Vi 'Vji which is Lorentz-invariant. Hence, symmetric choices 
are made and Qj.max = Qj.max = yJPi ■ Pj ■ Furthermore, the energy of the parton which 
initiates the shower is fixed to -Emax = Qmax^ in this way, one identifies the HERWIG 
showering frame. Setting < Q^^x ^^^t emission yields ( < 1 {9 < tt/2 in the 

massless approximatiom) in the showering frame. 

Although in the following we shall just deal with soft or collinear parton emissions, 
we point out that radiation of hard and large-angle partons can be implemented in 
HERWIG by applying matrix-element corrections. The region ( > 1, corresponding 
to hard and large-angle emissions, is the so-called 'dead zone' of the shower, where the 
exact amplitude is applied. Moreover, the 'hardest-so-far' emission in HERWIG showers 
is simulated according to the exact matrix element, thus allowing a smooth transition 
through the boundary of the dead zone [28]. 

Before discussing the implementation of medium effects, we wish to make a few 

^As discussed in [12], the z-limits in PYTHIA, for a branching with virtuahty p^, are given by 
Zmin — Po/p^ and Zmax = 1 " z^in, Pq bcing the infrared cutoff. This yields on average a larger 
2;-evolution range for PYTHIA. 

^Such a choice implies, e.g., that for e^e^ qq annihilation at energy ^/s, the initial value of the 
ordering variable for q and q is Qmax = y/s/2. 
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comments on the difference between HERWIG and PYTHIA showering algorithms. As 
discussed in the introduction, PYTHIA includes colour coherence only partially: in [12], 
medium modifications were included in virtuality-ordered showers, with an option to 
reject emission which do not fulfil angular ordering. As for matrix-element corrections, 
PYTHIA has no dead zones: it uses the soft/coUinear approximation throughout all 
physical phase space and simulates only the first branching via the exact amplitude 
[29]. Hadronization is finally implemented in PYTHIA following the string model [20]. 
Since parton showers and possible matrix-element matching are implemented in a pretty 
different way, if one studies parton-level quantities, as will be done hereafter, HERWIG 
results should not necessarily agree with the PYTHIA ones presented in [12], even 
though one models medium effects in the same way. 

Of course, we are aware that a comparison of the results yielded by HERWIG and 
PYTHIA is compulsory. However, it can be meaningfully done only at hadron-level 
and after both generators are tuned to the same data. In this work, we shall restrict 
ourselves to comment on possible discrepancies at parton level and leave to future work 
a more detailed comparison, which will demand the redoing of the previous extensive 
studies in hadron-hadron collisions. 



3 Medium-modified splitting functions 

In this section we discuss the main issues concerning the implementation of medium 
effects in parton shower algorithms. We follow the method presented in [10, 11] and add 
to the Altarelli-Parisi splitting function in the vacuum P{z) a medium-dependent term 
AP{z,p\E,q,Ly. 

P{z)^P{z) + AP{z,p\E,q,L). (3) 

In ([3]), is the branching-parton virtuality, E its energy, L the medium length, q the 
transport coefficient, defined as the average transverse momentum transferred from the 
medium to the parton per unity of free path [3]. Hereafter, we shall often make use of the 
quantities qL, the so-called accumulated transverse momentum, and of the frequency 
= qL^/2. 

The transformation ([3]) will be applied to the branching algorithm ([T]) and, in particu- 
lar, to the integrand function of the Sudakov form factor ([2]), taking care that HERWIG 
evolution variable is not the virtuality p^, but the energy- weighted angle defined 
above [fl. As in [12], a crucial hypothesis on which our work is based is that, even in 
a dense medium, the factorization ([1]) between branching and no-branching probability 
still holds and that the showering evolution variable is the same as in the vacuum. In 

•^For soft/coUinear radiation, it is straightforward to show that the relation ~ 2z(l — z)Q'^ holds. 
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HERWIG, this means that we shall assume angular ordering and colour coherence even 
for parton cascades in a medium. In other words, the colour flow will not be modified 
with respect to its vacuum pattern. 

Following [11], we shall modify the splitting functions for the branchings q qg 
and g gg, while we shall assume that the splitting g qq gets negligible medium 
modifications^ ^P{z,p'^, E,q, L) can be expressed in terms of the medium-induced 
parton radiation spectrum computed in the so-called BDMPS approximation 

[30-32]: 

APiz,p\E,q,L) ~ ^ , , (4) 

as dzap'^ 

As discussed in [12], Eq. (jlj) gives the medium corrections to the soft-divergent part of 
the Altarelli-Parisi splitting functions; the finite terms are assumed to be vacuum- 

AP{z,p'^, E,q, L) can be expressed in terms of the energy E of the parent parton 
and of the dimensionless variables E/uc and = k^/{qL), where /c^ is the transverse 
momentum of the radiated soft/collinear parton with respect to the splitting one. For 
multiple radiation, E will always be the energy of the splitting quark/gluon. As for 
the medium length, if Lq is the length for the first splitting, for the following emissions 
we have to take into account that a radiated parton, with energy fraction z, travels 
for a distance 2zE/k'^, the so-called parton formation length, before splitting again. 
Therefore, the effective medium length for subsequent branchings will be: 

L = Lo-^. (5) 

We note that in Eq. ([5]) L is not positive definite and, especially for small values of 
Lq and very soft/collinear splittings, it may well become negative after few emissions. 
Whenever this is the case, we shall assume that the shower continues in a vacuum-like 
fashion. 

In Ref. [11] the modified Altarelli-Parisi splitting functions were implemented in a 
Sudakov form factor computed in the same way as PYTHIA does. It was indeed found 
that medium effects have a remarkable impact and, for any given values of q and L, 
the Sudakov form factor is suppressed with respect to the vacuum one in the full p^- 
range, for both quarks and gluons. As a decreasing of the form factor corresponds to an 
enhancement of the branching probability, a higher parton multiplicity in the medium- 
modified shower must be expected. This was indeed found in [12] using Q-PYTHIA: 



^In any case, the splitting g — > qq is quite uncommon in the shower, since the spHtting function 
Pggi^) = Tfl[z^ + (1 — z)^] is not soft-enhanced. It is only after turning cluster hadronization on [19] 
that non-perturbative g qq splittings are forced and Pqg{z) plays a role. 

^That means that, e.g., for a, q ~> qg branching, where the Altarelli-Parisi splitting function is given 
by Pgq(z) — Cf[^ + (1 — medium-induced modifications affect the overall ^ 1/z factor, while 

the finite-z corrections remain unchanged. 
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although for a consistent comparison with data one must turn hadronization on, such 
an observation is in quahtative agreement with the expectations from radiative energy 
loss [1]. 

In Fig. [1] we present the HERWIG Sudakov form factor, possibly including medium 
effects by means of Eq. ([3]); for simplicity, we plot A(Q^, Q^ax) with fixed E and L = Lq, 
as if we were dealing with the the first emission in the shower. We consider medium- 
induced showers initiated by gluons of energy E = 10 and 100 GeV and media with q = 1 
and 10 GeV^/fm, Lq = 2 and 5 fm. In the following, for the sake of brevity, we shall 
label such medium configurations in terms of the accumulated transverse momentum, say 
qLo = 2, 5, 20 and 50 GeV^; the corresponding values of q and Lq are left understood. 
For gluon- initiated showers, the lower value of the evolution variable is, by default, 
Qo ~ 1.7 GeV [5], whereas we set by hand the upper valu^ Qmax = E. 

From Fig. [1] we learn that medium-induced effects are pretty relevant: for both 
E=10 and 100 GeV, the suppression due to including AP{z,p'^, E,q, L) can be up to 
several orders of magnitude, especially for relatively small values of Q. Of course, when 
Q approaches Qmax, all form factors tend to 1, and therefore the discrepancy due to 
the different values of qLo tends to become smaller, but it is nonetheless well visible 
throughout all Q-range. As for the behaviour of the modified As{Q'^, Qmax) when q and 
Lq change, we find that, for both E values and in all Q-range, the highest no-branching 

^The maximum value for the virtuality-ordered PYTHIA showers used in [12] is p^ax — 4i?^. Hence, 
even the p^-evohition range in PYTHIA is larger than the one in HERWIG. 




Figure 1: Gluon Sudakov form factors, in the vacuum (solid, black) and in media with 
qLo = 2 (dashes, blue), 5 (dots, red), 20 (dot-dashes, green) and 50 (magenta, solid) 
GeV^. The shown cases correspond to starting energies of = 10 GeV (a) and 100 
GeV (b). (For the interpretation of the references to colours in all figure legends, the 
reader is referred to the online version of this paper). 
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probability is obtained for qLo = 2 GeV^ and the lowest for qLo = 50 GeV^, as one 
should indeed expect. The comparison between the other two options in instead more 
cumbersome: for E = 10 GeV, Asi^Q"^, Qmax) for qL^ = 20 GeV^ is much larger than the 
form factor for (jLq = 5 GeV^, thus implying more branching in the case of accumulated 
transverse momentum equal to 5 GeV^. For E = 100 GeV, the gLg = 20 GeV^ form 
factor is above the qLo = 5 GeV^ one at low Q, but for Q values larger than 70-80 GeV, 
which roughly correspond to the first branchings in the shower, they are very close 
to each other|j Therefore, we can already foresee some similarities between these two 
medium parametrizations in the phenomenological analysis which we shall carry out for 
E = 100 GeV. In any case, the splitting-parton energy and the medium length vary 
throughout the cascade and Fig. [T] has been instead obtained for fixed Lq and E; the 
plotted As{Q^, Qmax)' though very useful to give us a qualitative estimate of the medium 
suppression, are not exactly the ones which we shall implement within the HERWIG 
algorithm. 

Since the ordering variables in PYTHIA and HERWIG are different, the numeri- 
cal comparison between Fig. [1] and the Sudakov form factors presented in [11] is not 
straightforward. However, as we said before, since the evolution range for variables z 
and in HERWIG is narrower than for the corresponding PYTHIA quantities, at 
least in the default scenarios, we can already predict a stronger medium-induced Su- 
dakov suppression and more branchings in PYTHIA rather than in HERWIG, in both 
vacuum and medium-modified showers. 



4 Results 

In this section we will present some phenomenological results showing the impact of 
medium modifications in HERWIG. In [12], study, the authors considered a 

gluon with fixed energy, giving rise to a medium-modified cascade within the PYTHIA 
model. This is not a standard option in HERWIG, but nevertheless, for the sake of 
comparison, we implemented a fictitious process with a single gluon of given energy E 
initiating a shower. In PYTHIA virtuality evolution, this was relatively straightforward, 
whereas in HERWIG, since one has angular ordering, in principle one would need two 
colour-connected partons to define the initial value of the evolution variable Q^ax- Still, 
we managed to start with a gluon of energy E, setting by hand the upper value of 
the evolution variable Q^ax = 2E'^. Subsequent splittings {g gg and so forth) will 
follow standard colour coherence 0. Similar tricks would be necessary if one wanted to 

^Note that the radiative energy loss is not only function of the product qLq, but its dependence on 
q and Lq is more complex. 

^As an alternative option, one can run the fictitious process e^e^ — > gg at ^/s = 2E, available in 
the HERWIG library, and study the parton shower only in one hemisphere. 
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simulate such a fictitious process in PYTHIA, but with showers ordered in transverse 
momentum, as their implementation follows a dipole formalism and one would need to 
start with two colour-connected partons [16]. Furthermore, as we are not investigating 
an actual physical process, we turn matrix-element corrections off; how to implement 
exact higher-order corrections to the hard-scattering process in a medium is anyway an 
open issue. 

As in [12], we shall study the distributions of the transverse momentum px of final- 
state partons, the emission angle 6, with respect to the initial-gluon axis, and the log- 
arithmic energy fraction, defined as ^ = ln(£'/|p|), \p\ being the modulus of the parton 
momentum. As we said before, we consider showers in the vacuum and in a dense 
medium characterized by parameters q and Lq. We shall assume gluon energies E = 10 
and 100 GeV and, as in the previous section, medium parameters leading to qLq = 2, 
5, 20 and 50 GeV^. Throughout the cascade, the effective medium length L will be ob- 
tained by applying Eq. ([5]); our distributions will be labelled in terms of the accumulated 
transverse momentum qLq. 

Though being an unphysical quantity, it is interesting to compute first the average 
parton multiplicity (A^) in the vacuum and in media characterized by the parameters 
given above. Such multiplicities strongly depend on the shower cutoff Qo and obviously 
increase whenever Qq is lowered. In Table [1] we quote these numbers for default HER- 
WIG and different values of qLo: as expected, parton multiplicities increase in a dense 
medium, which is in agreement with the measurements of larger medium-induced energy 
loss. For E = 10 GeV, the enhancement runs from 20% (gLo = 2 GeV^) up to about 
80% (qLo = 50 GeV^); for E = 100 GeV the corresponding numbers are 7% and 70%. It 
is interesting to notice in Table [1] that the cases qLo = 5 and 20 GeV^ give quite similar 
results and that the (jLq = 20 GeV^ multiplicity is slightly above the qLq = 5 GeV^ 
one for E = 10 GeV and below for a gluon of 100 GeV. In fact, the modified splitting 
functions (jl]) do not scale with respect to (jLq and do depend on the initiating parton 
energy: therefore, the phenomenological implications of the introduction of medium 
effects depend on the three variables q, Lq and E and different behaviours are to be 
expected if we start the cascade with partons of different energies, even though we use 
same values of q and Lq. In Fig. ([1]), although the Sudakov form factors were computed 
for fixed E and L, we already noticed differences between the options gLg = 5 GeV^ 
and (jLq = 20 GeV^, according to whether we have _E = 10 or 100 GeV. 

In Figs. |2H1] we present the pr, and ^ distributions, according to HERWIG in 
the vacuum and in a medium with the parameters q and Lq stated before, for energies 
E = 10 GeV (a) and 100 GeV (b). We plot everywhere normalized parton multiplicities, 
such as (1/A^) dN/dpT, so that the total integral will be equal to l|^ A feature of the 

^On the contrary, the plots presented in Ref. [12] arc normahzcd to the mean parton multiphcity; 
rescahng our distributions to the muhiphcities in Table [T] is straightforward. 



9 



Table 1: Average parton multiplicities in HERWIG showers initiated by gluons of energy 
of 10 and 100 GeV, in the vacuum and in a medium with assigned values of qL^. 



E 


qLo = 


qLo = 2 GeV^ 


qLo = 5 GeV2 


qLo = 20 GeV^ 


qLo = 50 GeV^ 


10 GeV 


2.56 


3.05 


4.14 


3.60 


4.56 


100 GeV 


6.95 


7.41 


8.79 


8.93 


11.70 



vacuum spectra at 10 GeV is that they exhibit a sharp peak at px = ^ = 6 = 0, 
corresponding to events where the initiating gluon evolves down to the shower cutoff 
with no branching at all. The probability for this to happen depends on the Sudakov 
form factor As{Ql^g^^,Ql), hence it may be eventually lowered by decreasing Qo or 
increasing Qmax, thus allowing more radiation in the evolution. Also, such a spike will 
vanish once hadronization is turned on. A peak, though less sharp, is visible even for 
media with small values of qLo, e.g. qLo = 2 GeV^, whereas it disappears for larger 
values of qLo- Also, it is less visible for E = 100 QeV, as in this case Qmax — 141.42 GeV 
and A((5max' Qo) a pretty small non-emission probability^"! It is interesting to notice 
that, even for ii^ = 10 GeV, such a peak was not so visible in [12], where the authors did 
find some events with no radiation, but, since the branching probability in the PYTHIA 
model is higher that in HERWIG in all evolution range, the Q-PYTHIA spectra exhibit 
smooth behaviour, even when p^, or approach zero. Later in this section, as an 
example, we shall plot integrated parton multiplicities, which behave smoothly also in 
the very first bin of our histograms. 

As for the px spectra, i.e. Fig. [21 for £'=10 GeV the vacuum distribution is above 
the others at large pt, as in [12]; however, due to the presence of the pt = peak, such a 
behaviour is less evident than in PYTHIA. In both cases, it is anyway the qLo = 5 GeV^ 
curve the one that gives the lowest normalized multiplicity at large px- At low pt, Q- 
PYTHIA gives the highest parton multiplicity in the two scenarios corresponding to 
L = 5 fm; the lowest multiplicity was instead the one for the vacuum case [12]. In 
HERWIG the low-p^ distribution is clearly affected by the px —>■ behaviour: the 
spectra exhibit a maximum about pr = 1 GeV, but, up to pt — 0.8 GeV, the difference 
between the four chosen qLo is quite small. For E = 100 GeV, the absence of a significant 
event fraction with no showering makes it easier to compare with the results yielded 
by Q-PYTHIA. Overall, for the reasons discussed above, namely lower effects in the 
Sudakov form factor, the impact of the implementation of medium modifications in 
HERWIG looks smaller than in Q-PYTHIA. However, the behaviour with respect to 

^"A similar effect was found, e.g., in Ref. [33], where the simulation of W/Z transverse momentum at 
the Tevatron exhibited a sharp peak ai px = 0, which disappears whenever one sets a non-zero intrinsic 
transverse momentum for the incoming partons. 
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1 2 3 4 5 10 20 30 40 

Pt (GeV) Pt (GeV) 

Figure 2: Transverse momentum multiplicity for a medium- modified shower initiated 
by a gluon of energy 10 GeV (a) and 100 GeV (b), in the vacuum (solid, black) and in 
media with accumulated transverse momentum qL^ = 2 (dashes, blue), 5 (dots, red), 
20 (dot-dashes, green) and 50 (solid, magenta) GeV^. 

varying q and Lq is qualitatively the same. The vacuum spectrum is the lowest at small 
Pt and the highest at large px, where the qLo = 50 GeV^ distribution is the one going 
to zero more rapidly. The other three curves, corresponding to qLo=2, 5 and 20 GeV^, 
lie in the range of the two extreme cases, as one should expect. In particular, at large 
Pt, qLo = 5 and 20 GeV^ seem to give similar distributions, whereas the (jLq = 2 GeV^ 
option yields the largest high-p^- normalized multiplicity after the vacuum. 

Concerning the 9 spectra, i.e. Fig. [3], first of all we need to point out an essen- 
tial difference between HERWIG and PYTHIA angular distributions, independently of 
medium effects. PYTHIA, which does not systematically implement angular ordering, 
allows parton radiation in all allowed phase space, i.e. up to = tt; on the contrary, 
in HERWIG angular-ordered showers, radiation is possible only for C < 1, ^ being the 
showering variable defined in Section 2. This leads to an empty region in the physical 
phase space, corresponding to large-angle emission. More precisely, in the massless ap- 
proximation, the condition ( < 1 implies 6 < tt/2 in the HERWIG showering frame. At 
the end of the showering, however, when jets acquire mass, a Lorentz boost is applied: 
this enlarges the angular permitted region, but a zone where radiation is forbidden is 
still present. In a physical process, such as e^e~ — ^ gg, this empty region would be filled 
by the soft radiation from the other jet, i.e. from q if one triggers the shower initiated 
by q 0. However, as we are simulating an unphysical process with only a single coloured 
parton, a dead zone, corresponding to large-angle radiation off the primary gluon, must 
be expected. 

^-'^In any case, even with two jets, there will still be a missing region, due to hard and wide-angle 
radiation, that one should eventually fill by means of matrix-element corrections. 
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In fact, the spectra in Fig. [3] show that the parton multiphcity above ~ 2 is 
neghgible. As medium modifications mostly affect the first few emissions, which, because 
of angular ordering, are the ones at the largest 9, it is reasonable that the spectra in a 
dense medium are above the vacuum ones in 1 < 6* < 2, as we learn from Fig. [31 At 
small 6', i.e. far from the hard scattering, the parton formation length 2zE/k^ is likely 
to be large and the effective medium length L, given by Eq. ([5]), becomes negative. This 
implies that small- angle branchings are typically vacuum-like. As far as the different qLo 
options are concerned, the qLo = 2 GeV^ is clearly the closest to the vacuum prediction, 
the difference between qLq = 5 and 20 GeV^ is small, though visible, whereas the 
qLq = 50 GeV^ spectrum, corresponding to the strongest medium modifications, leads 
to the highest parton multiplicity in the allowed large-angle range, say 1 < 6 < 1.8. 




e 



Figure 3: As in Fig. [21 but showing the angular distributions for showers initiated by a 
gluon of energy 10 GeV (a) and 100 GeV (b), in the vacuum and in a dense medium. 




1 2 3 4 0123456 



Figure 4: As in Figs. [2l and [3l but presenting the logarithmic energy-fraction (^). 
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We finally comment on the logarithmic energy-fraction {^) plots, presented in Fig. HI 
As for the E = 10 GeV case, the ^ = peak clearly spoils the low-^ part of the spectrum. 
However, regardless of such a peak, the suppression of partons with small energy fraction 
in a dense medium with respect to the vacuum is still clear. The behaviour of the four 
medium-modified spectra in terms of q and Lq is alike the one displayed in [12]. At 
small ^, the qLo = 2 GeV^ spectrum yields the highest energy fraction, followed, in 
order, by qLo = 20, 5 and 50 GeV^. Around the peak, the order is reversed, with 
qLo = 50 GeV^ giving the largest event fraction and qLo = 2 GeV^ the lowest. At 
large ^, the difference due to the different options of qLo tends to become smaller; at 
very high ^, the qLo=5 GeV^ spectrum is above the others. For E = 100 GeV, since 
the events without emissions are very few, the comparison among the different curves 
and the results in [12] is straightforward. The low-^ suppression in a medium with 
respect to the vacuum is still remarkable: within the modified spectra, at small ^, the 
qLo = 2 GeV^ option yields the highest event fraction, followed by qLo = 5 and 20 GeV^, 
whose results are roughly similar, and finally qLo = 50 GeV^. At large ^, the order of 
the different qLo cases is reversed, as happens for £"=10 GeV. 

The presence of a peak for a gluon with E = 10 GeV, whenever pt, or 6 are close 
to zero, can be a bit disturbing, since we have a non-negligible fraction of events in the 
very first histogram bin, affecting overall the full spectrum. As we said, a possibility 
to reduce the events in the first bin consists in enlarging the evolution range or turning 
hadronization on. As an alternative, for the time being, we can compute integrated 
quantities, such as, e.g., parton multiplicities up to a given value of ^, namely: 



In Fig. [5] we plot N{^) in the vacuum, for the four chosen scenarios of q and Lq, and 
E = 10 GeV: we observe a well visible medium-induced suppression at small ^ with 
respect to the vacuum. At low energy fractions, e.g. C, — 0.5, we find that the medium 
suppression is about 10% for qLo = 2 GeV^ and can be up to 60% for accumulated 
transverse momentum qLo = 50 GeV^. As we said in the introduction, although an 
actual comparison would require hadronization corrections and possible tuning of the 
Monte Carlo parameters to the data, such results are indeed quite encouraging, as they 
are in qualitative agreement with the small-^ suppression observed at RHIC. 

Before closing this section, following [12], we wish to study the pr, and ^ distribu- 
tions for fixed values of the medium length throughout the parton cascade. This means 
that, instead of applying Eq. (JSj) at each branching, we shall always employ L = Lq. In 
fact, because of Eq. ([5]), the medium length could even become negative at some point, 
which implies that, after the first few branchings, emissions occur in a vacuum-like fash- 
ion. For L = Lq, therefore, the medium length will always stay positive, and stronger 
medium effects are to be expected. As a case study, we investigate the two options 




(6) 
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12 3 4 

Figure 5: Integrated ^ distribution in the vacuum (solid, black) and in a medium with 
gLo= 2 (dashes, blue), 5 (dots, red), 20 (dot-dashes, green) and 50 (solid, magenta) 
GeV^. The energy of the gluon initiating the parton shower has been set to 10 GeV. 

qLo = 2 and qLq = 50 GeV^, for variable and fixed medium length, and still E = 10 and 
100 GeV. We first calculate the average parton multiplicity (N): we obtain, for L = Lq 
and E = 10 GeV, (A^) ^ 3.17 and 4.92, for qLo = 2 and 50 GeV^, respectively. This 
corresponds to enhancements of 4% and 8% with respect to the variable-length results 
in Table [H For an initiating gluon of energy 100 GeV, the average multiplicities are 
{N) ~ 8.09 {qLo = 2 GeV^) and {N) ~ 19.23 {qLo = 50 GeV^), hence the enhancement 
runs from 10% to 65%. Such results are in agreement with the expectation that a fixed 
length should emphasize medium-induced effects. 

A more evident impact of medium modifications is also exhibited by the pt, and 
^ normalized spectra, presented in Figs. EHHl In fact, for L = Lo, the suppression at 
large pt and the enhancement at small pt, displayed by Fig. [6l is stronger than what 
found when varying L according to Eq. ([5]). The effect is milder for qLo = 2 GeV^, 
as we are starting from small medium parameters, whereas it gets quite remarkable for 
qLo = 50 GeV^ and E = 100 GeV. In the angular distributions (Fig. [7j), for qLo = 
2 GeV^ the effect of fixing L is small but visible, while, for qL = 50 GeV^, the fixed- 
length angular spectrum is above the variable-L one for 0.4 < 6 < 1.2 and below at 
larger angles. In fact, setting L = Lq, quarks and gluons always 'see' a positive medium 
length, and partons, which in the variable-L cascade would be potentially produced at 
large 6, are now allowed to further branch. This explains why, for fixed L, we have a 
smaller normalized multiplicity at very large angles and a higher one at middle ^-values. 
As all distributions are normalized to unity, the fixed-L spectrum is above the variable-L 
one at small angles. Stronger suppression (enhancement) at small (large) values of ^ is 



14 




Figure 6: Transverse momentum distribution, with variable and fixed medium length 
L. The considered cases are qLq = 2 GeV^ (black, solid: variable L; dashes, blue: fixed 
L) and qLo = 50 GeV^ (red, dots: variable L; green, dot-dashes: fixed L), with an 
initiating gluon of 10 GeV (a) and 100 GeV (b). 




Figure 7: As in Fig. [6], but displaying the angular distributions for variable and fixed 
medium length. 



5 Conclusions 



We implemented medium-modified splitting functions in the HERWIG parton shower 
algorithm, whose evolution satisfies the angular ordering prescription. Following the 
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Figure 8: ^ spectrum for fixed and variable L. In the plots we follow the conventions 
adopted in Figs. [6] and [71 

Q-PYTHIA implementation detailed in [12], we added to the Altarelli-Parisi splitting 
function a term depending on the medium properties, the virtuality and the energy of the 
radiating parton. Such a modification was consistently implemented in both branching 
probability and Sudakov form factor. 

For the purpose of our phenomenological analysis, we included in HERWIG a ficti- 
tious process, where a shower is originated by a single gluon with energy 10 and 100 
GeV. We studied the usual vacuum case, as well as showers in a medium characterized 
by transport coefficient q and length L, with L varying throughout the parton cascade. 
As HERWIG showers satisfy angular ordering, when starting with just one parton and 
not with a pair of colour-connected ones, we had to set by hand the upper value of the 
evolution variable. 

We ran HERWIG with modified splitting functions and, first of all, observed that 
the average parton multiplicity in a medium is higher than in the vacuum, with an en- 
hancement which varies from 20% to 80% when the accumulated transverse momentum 
qLo runs between 2 and 50 GeV^. Then, we studied differential distributions, such as 
transverse momentum (pt), angle [6) and logarithmic energy- fraction (^) spectra. In- 
deed, we found the features which one should expect in a dense medium with respect to 
the vacuum: enhancement at small px and suppression at high pt, more partons at large 
angles, suppression at low ^, compensating an enhancement for middle-high ^ values. 
As for the results obtained for E = 10 GeV, we noticed a sharp peak aX px = = ^ = 0, 
corresponding to a fraction of events with no emission at all. For this reason, we also 
plotted the ^ integrated distribution, which is insensitive to such a peak, and found re- 
sults more similar to the ones yielded the Q-PYTHIA code and in qualitative agreement 
with the observations at RHIC. We investigated the differential spectra for fixed values 
of the length and found stronger medium effects, as in this case L does not decrease 
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throughout the shower and is always positive. 

As for the comparison with the results in [12], overall we found acceptable quali- 
tative agreement, although, due to the fact that the two shower algorithms are pretty 
different and that we have been considering unphysical parton-level quantities, with 
hadronization switched off, some quantitative discrepancies are well visible. In fact, 
we observed that due to the choice of the showering variables, more branchings are 
simulated in PYTHIA, at least when using the default parametrizations, even before 
including medium modifications. In particular, sharp zero-peaks, exhibited by all HER- 
WIG distributions when considering a starting gluon of 10 GeV, were not present in 
PYTHIA. 

We believe that it will be now very interesting to consider an actual physical process, 
e.g. at the LHC, investigate hadron-level observables and eventually understand how, 
after including medium sphtting functions, modified HERWIG fares with respect to Q- 
PYTHIA. For such a comparison to be trustworthy, both generators should be tuned 
to the same data set. Also, it will be very useful understanding whether the medium- 
induced effects found in this analysis still persist once hadronization is turned on. As 
for the issue of angular ordering and evolution variables, it will be cumbersome studying 
quantities sensitive to colour coherence or non-global observables, as in Refs. [15,17], and 
understand whether the discrepancies between HERWIG and PYTHIA still persist or 
medium-effect implementations wash them out. However, as we said in the introduction, 
a comparison with hadron-hadron and heavy-ion data is compulsory before making any 
strong statement on the role played by the shower evolution variable in a dense medium. 
The documentation of the Q-HERWIG code, including medium-modified splitting func- 
tions, which should be not seen as an official release, but rather as an add-up to the 
latest fortran version [5] , is currently in progress [24] . 
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